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1.  SUMMARY  AND  INTRODUCTION 

1 . 1 SUMMARY 

A study  was  performed  to  evaluate  failure  probabilities  for  the 
cylindrical  tunnels  in  the  DIABLO  HAWK  test,  to  recommend  test  changes, 
and  to  investigate  the  effects  of  correlation  on  the  failure  probabilities. 
It  is  anticipated  that  either  5 or  6 tunnels  will  fail  in  the  DIABLO  HAWK 
test,  with  5 tunnels  experiencing  severe  damage.  A brief  decision  analysis 
indicates  that  the  test  would  be  optimized  with  the  addition  of  one  struc- 
ture designed  and  situated  in  order  to  have  a failure  probability  of  .5. 
However,  the  new  economic  benefit  is  small  and  in  the  long  run,  the 
omission  of  the  additional  structure  will  have  very  minor  influence.  If 
the  structural  properties  are  fully-correlated,  and  if  the  applied  loading 
(on  the  several  tunnels)  is  also  fully-correlated,  then  the  predicted 
number  of  failures  for  DIABLO  HAWK  is  expected  to  fall  between  3 and  8 
with  90%  probability.  Additional  results  of  the  study  include  the  con- 
clusion that  fragility  curves  can  readily  be  developed  for  hollow  spheres, 
and  the  elasto-plastic  design  of  spheres  is  amenable  to  computer  implemen- 
tation but  is  somewhat  more  complicated  than  tunnel  design. 

1.2  INTRODUCTION 

Previous  studies  (Reference  1)  were  reasonably  successful  in  pre- 
dicting the  number  of  tunnel  failures  in  the  MIGHTY  EPIC  test.  Never- 
theless, the  previous  analyses  did  not  address  the  DIABLO  HAWK  test,  nor 
did  they  include  the  important  effect  of  correlation  among  the  random 
variables  involved. 

The  designers  of  the  (MIGHTY  EPTC/DIABLO  HAWK)  test  structures 
wanted  to  investigate  the  effects  of  multiple  loadings  and  also  the 
propagation  of  a ground  shoc)c  pulse  through  ”pre-shoc)ced"  material.  Neither 
of  these  effects  are  addressed  in  the  present  report.  If  one  neglects 
these  effects  (as  the  present  report  does)  then  it  appears  that  the  DIABLO 
HAWK  test  is  less-than-optimally  designed.  Suggestions  are  made  herein 
to  improve  the  probability  of  obtaining  useful  data,  without  altering  the 
effects  of  multiple  loading  or  pre-shoclcing . 

The  two  prime  sources  of  correlation  for  the  tunnel  structures  are 

• The  commonality  of  construction  procedures,  concrete  strengths, 
etc.  that  may  malce  the  structural  properties  highly  correlated, 
and 

3 


J 


• The  correlation  of  weapon  output  from  structure  to  structure, 
namely,  if  the  free-field  pressure  is  high,  it  will  be  high 
for  all  the  structures,  etc. 

The  simplest  case  of  full  correlation  is  the  easiest  to  handle  computa- 
tionally. But,  in  fact,  since  the  rock  properties  vary  throughout  the  test 
site,  the  "fully-correlated"  approximation  does  not  represent  the  real, 
physical  situation.  Thus,  it  is  expected  that  the  actual  test  will  lie 
somewhere  between  the  totally  uncorrelated  and  fully  correlated  extremes. 

From  physical  arguments,  one  can  determine  that  the  effect  of  correla- 
ting the  structural  properties  and  the  pressure  loading  will  be  to  emphasize 
the  extremes  of  the  probability  distributions.  This  intuitive  result  is 
borne  out  by  the  calculations. 

The  question  of  developing  fragility  curves  for  hollow  spheres  is 
also  addressed,  as  is  the  problem  of  elasto-plastic  design  of  spheres 
in  rock.  Both  of  these  problems  are  conceptually  straightforward,  although 
the  elasto-plastic  design  is  more  complicated  for  spheres  than  it  is  for 
cylinders. 

The  report  includes  major  sections  on  predictions  and  recommendations 
for  DIABLO  HAWK,  the  effect  of  correlation,  sphere  calcualtions,  and 
concluding  remarks. 
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2. 


PROBABILITY  OF  FAILURE 
FOR  DIABLO  HAWK  CYLINDERS 


2 . 1 BACKGROUND 

In  the  previous  report  (Reference  1)  probabilistic  predictions  were 
made  regarding  the  expected  number  of  failures  (for  lined  tunnels  in  rock) 
in  the  MIGHTY  EPIC  test.  Assuming  a factor  of  f 1.1  uncertainty  in  the 
free-field  pressure  produced  in  the  tuff,  it  was  stated  that  ‘the  number  of 
expected  failures  lies  between  4 and  6,  within  90  per  cent  confidence 
limits.”  In  the  test,  four  (4)  tunnels  were  damaged  significantly,  with 
two  of  the  four  suffering  severe  damage  and  two  suffering  moderate  damage. 

The  two  structures  that  were  severely  damaged  had  a mean  probability 
of  failure  P^  equal  to  .99,  calculated  for  the  factor  of  1.1  uncertainty. 
The  other  two  structures,  which  received  moderate  damage,  had  failure 
probabilities  P^  = .91  and  P^  » .79,  respectively.  It  is  noteworthy  that 
one  other  structure,  which  experienced  only  slight  damage,  also  had  a 
failure  probability  P^  = .91. 

In  the  paragraphs  which  follow,  the  calculation  of  similar  failure 
probabilities  for  the  DIABLO  HAWK  test  are  outlined.  The  results  show 
that  five  (5)  tunnels  have  failure  probabilities  P^  > .97  for  a factor  of 

1.1  uncertainty  in  the  free-field  stress.  Consequently,  it  is  predicted 
that  either  5 or  6 tunnels  will  fail  (within  90  per  cent  confidence  limits) 
in  the  DIABLO  HAWK  event.  Furthermore,  these  5 tunnels  are  expected  to 

be  sever ly  damaged,  in  a manner  like  the  two  MIGHTY  EPIC  structures  which 
had  failure  probabilities  > .99. 

The  reader  is  again  cautioned  that  these  predictions  do  not  include 
the  effect  of  a second  loading  on  the  rock  medium.  If  the  second  loading 
causes  major  effects  not  anticipated  at  this  time,  then  the  predicted 
number  of  failures  might  change  significantly. 

2.2  EFFECTS  OF  UNCERTAINTY  IN  FREE-FIELD  PRESSURE 

In  order  to  understand  the  effect  that  uncertainties  in  the  free- 
field  stress  (pressure)  have  on  the  mean  probability  of  failure  (P^) , it 
is  worthwhile  to  consider  the  steps  in  the  calculation  of  Pf  . This  pro- 
cedure is  outlined  in  the  flow-chart  of  Figure  2-1.  The  individual  blocks 
in  the  flow-chart  can  be  briefly  summarized  as  follows; 

Step  1:  Calculate  the  design  stress,  o^,  for  a particular  tunnel. 
The  design  stress  depends  upon  the  steel  thickness,  h^,  the  concrete  thick- 
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READ  DESIGN  RANGE.  Rq,  FROM 
PRESSURE  vs.  RANGE 
CHARACTERISTIC 


KNOWING  actual  RANGE.  R, 
COMPUTE  RATIO  R,„ 


ENTER  CURVE  FOR  P^  vs,  R^j, 
AND  READ  OFF  FAILURE  ° 
PROBABILITY 


STRUCTURES 


Figure  2-1.  Flow  of  The  Calculations 


ness,  t , etc.  A failure  strain  of  c.  - .05  was  used  (cf.  Reference  1) 
c 0 

along  with  mean  values  of  the  design  variables  E , v , etc.  The 

rock  rock 

mean  values  of  these  variables  are  listed  in  Table  2-1.* 

Step  2:  Using  the  design  stress,  o , and  a curve  of  free-field 

o 

pressure  vs.  range  (estimated  for  the  event) , one  can  read  off  a corres- 
ponding design  range,  R^. 

Step  3:  From  a knowledge  of  the  actual  range,  R (from  the  structure 

to  the  weapon  point)  and  the  design  langa,  R , from  Step  2,  one  can  compute 

the  dimensionless  ratio  (R/R  ) . 

o 

Step  4:  Using  the  curves  developed  previously,  giving  vs. 

(R/R^) , one  can  read  off  the  corresponding  mean  probability  of  failure 
(see  Figure  2-2) . 

Step  5:  This  procedure  is  repeated  for  all  the  structures  involved. 

The  calculations  just  outlined  were  p>erformed  for  the  DIABLO  HAWK 
test  cylinders  using  uncertainty  factors  of  2 and  4 for  the  free-field 
pressure.  Similar  (but  somewhat  simpler)  calculations  were  done  using 
an  uncertainty  factor  1.1.  The  results  of  these  calculations  are  given 
in  Table  2-2.  As  can  be  seen  by  referring  to  the  Table,  the  effect  of 
decreasing  the  uncertainty  in  free-field  pressure  is  to  decrease  the  prob- 
abilities when  P^  ^ .5  and  to  increase  the  probabilities  when  P^  ^ .5. 

This  sharpening  of  the  distribution  is  evident  in  Figure  2-3,  which  also 
shows  a suggested  "optimum"  test  design.  The  reader  should  be  aware 
that  the  curves  of  Figures  2-2  and  2-3  are  based  upon  an  assumed  value 
of  the  "failure  strain,"  namely  Eg  = .05.  Although  this  value  (Eg  * .05) 
gave  reasonably  good  results  for  the  MIGHTY  EPIC  test,  the  specification 
of  a failure  strain  is  an  involved  subject  and  deserves  further  study. 

The  curves  of  Figure  2-3  substantiate  the  previous  finaings  (Ref- 
erence 1,  p.  5-6)  that  "assuming  a factor  of  4 uncertainty  in  the  free- 
field  pressure,  the  actual  test  was  reasonably  well-designed."  Conversely, 
when  the  uncertainty  is  reduced  to  a factor  of  1.1  (which  is  thought  to 
be  representative  of  the  well-calibrated  Nevada  Test  Site)  it  appears  that 
the  combined  (MIGHTY  EPIC/DIABLO  HAWK)  test  cylinders  were  less  than 
optimally  designed.  This  finding  is  reinforced  when  the  probability 
of  exactly  N failures  is  comptuted,  as  discussed  in  the  following  paragraph. 

♦These  tabulated  values  were  used  previously  in  Reference  1. 
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Figure  2-2.  P.  - Maps  Using  Two  Different  Uncertainty  Factors  (From  Reference  1 ) 


Table  2-2  describes  11  individual  "events"  (in  the  terminology  of 
probability  theory)  each  with  a corresponding  failure  probability.  (Note 
that  at  this  point  we  are  treating  each  test  cylinder  as  independent  - 
i.e.,  with  no  correlation  among  the  structures  or  the  loading.  The  influ- 
ence of  correlation  is  discussed  in  Section  4.)  Viewed  as  11  individual 
events,  where  the  i—  structure  has  a probability  of  failure  P^,  Collins 
(Reference  2)  has  shown  how  to  compute  the  probability  of  exactly  zero 
failures,  etc.  The  mathematical  expressions  are  presented  in  References  1 
and  2,  and  calcualtions  are  most  amenable  to  a simulation  technique.  The 
results  of  this  simulation  are  shown  in  Figure  2-4  for  the  three  uncertainty 
factors,  namely  4,  2,  and  1.1. 

Note  that  Figure  2-4  (c)  for  a factor  of  4 uncertainty  shows  a 
symmetrical  histogram,  (indicating  a well-balanced  design)  whereas  Figure 
2-4 (a)  shows  that  (within  90  per  cent  confidence  limits)*  5 or  6 cylindrical 
structures  will  fail  in  the  DIABLO  HAWK  test.*  In  fact,  the  probability 
that  exactly  5 cylinders  fail  is  approximately  .65,  which  reflects  the  fact 
that  5 structures  have  individual  failure  probabilities  ^ .97.  These 
results  are  presented  in  Table  2-3  as  "Predictions  for  DIABLO  HAWK". 

As  indicated  by  the  footnote  in  Table  2-3,  thus  far  the  effects  of 
initial  ovaling  (i.e.,  initial  imperfections)  caused  by  the  MIGHTY  EPIC 
test  have  not  been  included,  nor  have  the  effects  of  correlation  among 
structures  or  loading  been  examined.  These  subjects  are  discussed  in  the 
sections  which  follow. 

2.3  EFFECTS  OF  INITIAL  OVALING  FROM  MIGHTY  EPIC  TESTS 

To  understand  the  effect  of  initial  ovaling  (initial  imperfections) 
on  the  behavior  of  lined  tunnels  in  roc)c,  it  is  instructive  to  consider  an 
analogous  problem,  namely  buc)cling  of  an  imperfect  column.  The  problem 
of  a slightly  croo)ced  column  is  discussed  in  Reference  3,  p.  13,  and  also 
in  Timoshenlco’ s well-)cnown  text  (Reference  4). 

For  a perfect  simply-supported  column,  of  length  L and  having  a 
uniform  bending  stiffness  El,  the  elastic  buc)cling  load  is  given  by 

^The  term  "confidence  limits"  is  discussed  in  most  texts  on  probability 
and  statistics  (e.g.  Reference  1) . 

♦This  statement  assumes  that  an  uncertainty  factor  of  1.1  is  repres'nta- 
tive  of  the  Nevada  Test  Site. 

♦♦These  results  are  adapted  from  Brush  and  Almroth,  Reference  3. 
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Figure  2-4.  Histograms  for  Diablo  Hawk  Test  (uncorrelated) 
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Table  2-3.  Predictions  for  Diablo  Hawk* 


STRUCTURE 

FAILURE 

NUMBER 

PROBABILITY 

^1.1. 

CX  10 

.00 

CX  9 

.00 

NOT  EXPECTED  TO  FAIL 

CX  8 

.00  . 

1 

CX  7 

.11 

POSSIBLE  ONE  OF  THESE  THREE  WILL 

CY  23 

.15 

FAIL.  WITH  DAMAGE  LIKE  SLIGHT  FLAT- 

CY  22 

.11  . 

TENING  EXPERIENCED  IN  MIGHTY  EPIC. 

CY  16 

.97 

j 

CY  15 

.98 

THESE  FIVE  ARE  EXPECTED  TO  FAIL. 

1 

1 

CY  14 

.9995 

EXPECT  SEVERE  DAMAGE  COMPARABLE  TO 

CZ  3 

.9995 

FAILURES  IN  MIGHTY  EPIC. 

j 

CZ  1 

.9995  . 

‘NOTE  THAT  THESE  PREDICTIONS  DO  NOT  INCLUDE  THE  EFFECTS  OF  IMPERFECTIONS, 
(E.G.,  INITIAL  OVALING  CAUSED  BY  MIGHTY  EPIC  TEST)  NOR  DO  THEY  INCLUDE 
THE  EFFECTS  OF  CORRELATION  AMONG  STRUCTURES  OR  LOADS 


For  a perfect  simply-supported  column,  of  length  L and  having  a 
uniform  bending  stiffness  El,  the  elastic  buckling  load  is  given  by 


When  the  column  is  initially  imperfect,**  say  it  has  an  initial  deflection 
w*  = sin  — (2-2) 

then  the  deflection  in  the  buckling  mode  is  given  by 
PC  * sin  ^ 

w(x)  = — 5 — (2-3) 

ILIEI  _ p 

where  Cj^*  is  the  amplitude  of  the  imperfection.  Note  that  equation  (2-3) 

predicts  that  the  lateral  deflection,  w,  at  the  middle  of  the  column  (i.e., 

2 2 

X = L/2)  becomes  unbounded  as  the  load  P approaches  tt  EI/L  . Qualitatively, 
the  effect  of  initial  imperfections  is  shown  in  Figure  2-5. 

When  the  initial  imperfection,  is  small,  the  behavior  approaches 

that  of  a perfect  column.  As  Cj^*  becomes  larger,  the  lateral  deflection 

grows  more  rapidly,  but  in  all  cases,  the  lateral  deflection  increases  with- 

2 2 

out  bound  as  P approaches  P^^  = tt  EI/L  , the  buckling  load  of  the  perfect 
column.  This  result,  namely  that  the  maximum  load  which  the  column  will 
carry  is  not  strongly-dependent  on  the  initial  imperfection  (C^*)  has  led 
researchers  to  classify  columns  as  "insensitive"  to  initial  imperfections. 

As  might  be  imagined,  not  all  structures  are  insensitive  to  initial 
imperfections.  Figure  2-6  qualitatively  shows  the  effect  of  imperfections 
on  the  buckling  of  an  axially  compressed  cylindrical  shell  (Reference  3 ) . 
Certain  shell  structures  are  notorious  for  their  sensitivity  to  initial 
imperfections,  and  this  fact  usually  is  manifested  by  a significant  dis- 
crepancy between  (perfect)  theory  and  (imperfect)  experiments.  With  respect 
to  the  problem  at  hand,  namely  failure  of  lined  tunnels  in  rock,  the  next 
step  is  to  determine  whether  the  tunnels  are  relatively  insensitive  (like 
the  column)  or  highly-sensitive  (like  the  axially-loaded  shell) . The  results 
available  (Reference  3 ) indicate  that  tunnels  loaded  by  external  (hydro- 
static) pressure  are  relatively  insensitive  to  initial  imperfections,  at 
least  for  the  (MIGHTY  EPIC/DIABLO  HAWK)  geometries. 
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Figure  2-6.  Load-Deflection  Curves  for  an  Axially-Compressed  Cylindrical  Shell 


For  example,  Figure  2-7  shows  the  comparison  between  theory  (the 
solid  line)  and  experiment  (the  open  dots)  for  buc(tling  pressure  p vs.  the 
Batdorf  parameter,  Z,  When  Z is  computed  for  the  test  cylinders,  one  has 


Z 


(1 


v2,1/2 


(2-4) 


where 

L = 16  ft  (length) 
a = 2 ft  (radius) 
h = .75  in  (thic)tness) 

V = .3  (Poisson's  ratio) 

which  gives 

Z 73  = 1953 

when  a steel  thiclcness  of  .75  inches  is  used. 

This  calculation  does  not  reflect  the  fact  that  the  steel  liner  is 
attached  to  the  concrete,  which  is  much  thic)cer.  By  equating  the  bending 
stiffness  (El)  of  the  steel-and-concrete  to  that  of  a completely  steel 
"equivalent"  shell,  one  can  conservatively  estimate  an  equivalent  thic)cness, 
hgg,  as  less  than  12  inches.  In  this  case,  one  computes 

Zj^2  = 122  (2-5b) 


Referring  now  to  Figure  2-7,  one  sees  that  the  Z values  of  equations 
(2-5)  fall  in  a region  where  theory  and  experiment  are  in  reasonable  agree- 
ment. This  agreement  is  usually  lac)cing  if  the  structure  is  sensitive  to 
initial  imperfections.  To  further  substantiate  this  point.  Figure  2-8 
shows  a theoretical  "imperfection  sensitivity  parameter",  a2,  plotted  vs.  Z 
for  cylindrical  shells  loaded  by  external  pressure.  Note  that  for  Z > 100, 
the  sensitivity  parameter  a^  is  small,  (tending  toward  zero)  which  indicates 
a lac)c  of  sensitivity  to  initial  imperfections. 

The  net  result  of  this  discussion  is  that  the  (MIGHTY  EPIC/DIABLO 
HAWK)  test  cylinders  are  not  expected  to  be  sensitive  to  initial  imper- 
fections. Thus,  although  MIGHTY  EPIC  produced  some  (relatively  slight) 
initial  ovaling  in  the  (undamaged)  test  cylinders,  these  imperfections 
are  not  epxected  to  significantly  change  the  DIABLO  HAWK  predictions. 

It  is  )cnown  that  the  effect  of  the  imperfections  will  be  to  increase  the 
expected  probability  of  failure,  P^,  but  such  changes  are  not  expected 
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7.  Comparison  of  Theoretical  and  Experimental  Values  for  Cylinders 
Subjected  to  Hydrostatic  Pressure 


Figure  2-8.  Imperfection  Sensitivity  of  Cylinders  Subjected  to  Lateral  Pressure 


to  be  significant.  Again  the  reader  is  cautioned  that  the  effect  of 
pre-stress  in  the  rock  medium  has  not  been  inlcluded  in  this  discussion 
and  this  pre-stress  effect  may  be  important. 


Referring  back  to  Figure  2-4  for  a moment,  one  sees  that  5 or  6 
cylindrical  structures  are  expected  to  fail  in  the  DIABLO  HAWK  test,  and 
Table  2-3  indicates  that  5 cylinders  will  be  "severely"  damaged.  Armed 
with  this  knowledge,  one  might  pose  the  question,  "How  can  we  adjust  the 
DIABLO  HAWK  test  structures  to  improve  the  probability  of  obtaining 
successful  test  data?"  This  question  is  addressed  in  Section  3.  of  this 
report. 
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3. 


RECOMMENDATIONS  FOR  THE 
DIABLO  HAWK  TEST 


3.1  RECOMMENDATIONS 

The  test  planners  for  the  (MIGHTY  EPIC/DIABLO  HAWK)  structures  test 
had  certain  design  goals  in  mind  when  they  selected  the  sizes  and  placement 
of  their  structures.  Although  all  of  their  design  goals  are  not  jtnown  to 
the  writer,  two  goals  which  were  desired  are 

• The  desire  to  load  the  same  structures  more  than  once,  i.e., 
multiple  loadings,  and 

• The  desire  to  propagate  a second  shoe)?  pulse  through  the  tuff/ 
grout  medium,  which  may  have  been  fractured  and/or  compressed 
by  the  first  test. 

These  two  test  goals  have  an  impact  on  the  realistic  (i.e.,  practical) 
changes  which  can  be  considered  to  improve  the  probability  of  acquiring  use- 
ful data  in  the  DIABLO  HAWK  test.  Section  2 tells  the  story  of  what  is 
expected  on  DIABLO  HAWK: 

• Five  cylinders  are  expected  to  fail  severely,  similar  to  those 
most  extensively  damaged  in  MIGHTY  EPIC,  and 

• One  other  cylinder  (out  of  three  candidates)  may  experience  slight 
buc)cling. 

Faced  with  these  goals  and  predictions,  the  following  recommenda' ion 
is  made: 

• Add  one  or  more  cylindrical  structures  with  a .5  probability  of 
failure. 

The  reasons  behind  this  recommendation  are  given  in  the  paragraphs  which 
follow,  along  with  other  alternatives  that  appear  less  practical.  Section 

3.2  also  contains  an  example  decision  analysis  which  also  indicates  the  need 
for  one  more  structure. 

Referring  to  either  Figure  2-3  or  Table  2-2,  for  a factor  of  1.1 
uncertainty,  one  sees  a less  than  optimally  designed  experiment  for  the 
cylindrical  tunnel  structures.  Five  structures  have  ^ .97;  three 
structures  have  P^  .00,  and  the  other  structures  have  mean  failure  prob- 
abilities of  .11,  .11,  and  .15  respectively.  To  improve  the  test  design, 
it  is  recommended  that  DNA  consider  adding  at  least  one  more  tunnel,  designed 
with  an  intermediate  failure  probability. 
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Such  a test  matrix  will  give  a more  optimum  design  and  is  expected 
to  increase  the  probability  of  obtaining  useful  data.  If  no  changes  are 
made  in  the  DIABLO  HAWK  test,  it  is  very  probable  that  five  (5)  structures 
will  be  severely  damaged,  one  just  slightly  damaged,  and  the  other  five 
will  be  basically  unaffected  by  the  event.  The  proposed  additional  structure 
is  expected  to  produce  an  intermediate  damage  which  will  increase  the  infor- 
mation and  useful  data  obtained. 
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other  alternatives  which  might  be  considered  but  which  have  inherent 
difficulties  are  listed  in  Table  3-1,  along  with  a qualitative  estimate  of 
their  attendant  cost.  It  is  evident  that  further  study  of  the  problem  is  i 

required  before  an  intelligent  final  choice  can  be  made  among  the  various  i 

alternatives.  In  that  regard,  it  may  be  premature  to  suggest  adding  more 
tunnels  to  improve  the  test,  but  the  additional  tunnels  appear  to  be  the 
leading  candidate  for  improvement  as  of  this  writing. 

The  reader  should  note  that  the  predictions  made  in  Section  2 and 
re-stated  herein  do  not  include  the  effects  of  "pre-shocking"  on  the  tufr. 

The  MIGHTY  EPIC  test  fractured  some  of  the  rock,  compressed  it,  etc,,  and  the 
predictions  of  this  report  have  not  included  these  effects.  It  is  assumed 
tnat  these  effects  are  negligible  or  of  minor  significance;  however,  this 
assumption  may  prove  to  be  incorrect. 

3.2  EXAMPLE  DECISION  ANALYSIS  FOR  SELECTING  AN  OPTIMUM  NUMBER  OF 

ADDITIONAL  STRUCTURES 

Although  the  DIABLO  HAWK  test  is  already  set,  and  further  changes  will 
be  difficult,  it  is  of  interest  to  determine  whether  additional  cylinders 
would  be  of  value  in  the  future  application  of  data  from  the  test.  The  problem 
with  the  current  configuration  for  DIABLO  HAWK  is  that  the  cylinders  are  either 
failing  with  a high  probability  or  with  a very  low  probability,  and  there 
is  very  limited  data  as  to  the  average  probability,  i.e.,  the  level  of 
structure  strength  where  the  probability  of  failure  would  be  around  .5. 

Consider  the  problem  shown  in  Figure  3-1.  There  is  a design  level  in  a 
future  system  which  resulted  from  the  conclusions  about  cylinder  strength  in 
the  MIGHTY  EPIC  and  the  DIABLO  HAWK  tests.  There  is  an  uncertainty  about 
this  mean  strength,  Oq,  resulting  from  the  lack  of  information  from  failures 
in  the  tests  relating  to  the  P^  = . 5 range.  Consequently,  a normal  distri- 
bution is  shown,  representing  the  uncertainty  in  the  true  capability  of  the  ; 

structure.  Since  this  uncertainty  exists  at  the  time  of  the  design  of  the 
future  system,  the  designer  must  take  a conservative  position  and  possibly 
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overdesign,  because  he  is  not  sure  of  the  true  mean  capability  of  the  cyl- 
inders. The  resulting  cost  to  the  future  system  is  in  additional  materials 
which  would  be  needed  to  give  the  additional  strength.  If,  indeed,  the 
cylinder  were  stronger  than  estimated  by  the  designer  because  of  his  lack 
of  uncertainty,  then  there  would  be  an  economic  opportunity  loss  (a  spending 
of  money  which  did  not  need  to  be  spent) , due  to  the  uncertainty. 


COST  OF 
MATERIALS 
IN  DESIGNED 
STRUCTURE 


DOLLAR  SAVINGS  IN  MATERIALS  IF 
DESIGN  COULD  BE  REDUCED 

.0  , UNCERTAINTY  IN  ACTUAL  CAPABILITY 
° RESULTING  FROM  UNCERTAINTY  IN 


PAST  TEST  RESULTS 


(SOFTER) 


1 .0 

-DESIGN  LEVEL 


(HARDER) 


■DESIGN  LEVEL  RESULTING 
FROM  CONCLUSIONS  REGARDING 
TEST,  u (IN  THIS  CASE  u = 1 ) 


Figure  3-1.  Economic  Opportunity  Loss  as  a Function 
of  Mean  Design  Level  and  Uncertainty 


The  average  financial  loss  due  to  the  uncertainty  can  be  expressed 
as  an  "expected  opportunity  loss"  as  shown  in  Equation  (3-1) 


E [Opportunity  Loss]  = EOL = 


/(x-u^) 

\2  0 2 
' o 


(3-1) 


Note  there  is  no  opportunity  loss  if  the  design  level  should  have  actually 
been  higher.  On  the  other  hand  this  would  increase  p£  which  is  also  undes- 
irable. Hence  another  line  should  start  at  x = and  have  a positive  slope 
based  on  some  financial  interpretation  of  increased  probability  of  failure. 
At  this  point  we  will  suggest,  for  purposes  of  demonstration,  that  this 
second  line  be  ignored  and  that  the  problem  be  restricted  to  that  shown  in 
Figure  3-1. 

The  uncertainty  in  the  proper  design  level  can  be  reduced  by  adding 
structures  to  the  next  test  which  are  designed  such  that  they  have  a prob- 
ability of  failure  equal  to  approximately  .5.  It  would  be  best  to  have  a 
range  of  failure  probabilities  from  ,2  or  .3  to  .7  or  .8,  but  regardless 
the  results  should  yield  a mean  design  level,  x,  which  would  represent  the 
best  estimate  of  the  average  design  level. 
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Assume  that  the  uncertainty  in  the  failure  level  of  any  specific 
structure  (Cg)  can  be  obtained  from  the  fragility  curve.  Then  the  standard 
error  of  the  estimate  of  the  mean  is  = o^/Zn,  where  n is  the  number  of 
structures  being  tested. 

To  determine  the  optimum  number  of  tests  required  to  minimize  cost 
(combined  system  cost  and  testing  cost) , assume  that  the  normal  distribution 
in  Figure  3-1  represents  our  "prior"  knowledge  of  the  uncertainty  and  theor- 
etically we  could  afford  to  pay  up  to  the  number  of  dollars  in  the  EOL  for 
additional  tests  ^ the  results  from  the  additional  tests  were  able  to  reduce 
to  zero  our  uncertainty  in  the  estimate  of  the  true  mean.  Unfortunately,  the 
results  of  the  tests  will  still  produce  uncertainty  in  the  estimate  of  the 
mean,  although  less  than  the  original. 

If  it  is  assumed  that  all  the  distributions  involved  are  normal  (i.e.  , 
Gaussian) , then  it  can  be  shown  that  the  revised  estimate  of  the  uncertainty 
in  the  estimate  of  the  true  mean  can  be  written  as  follows  (Ref.  14,  15) 


where  a is  the  revised  uncertainty  in  the 
^ estimate  of  the  true  mean 

is  the  prior  uncertainty 

o-  is  the  uncertainty  in  x resulting 
^ from  sampling  (additional  tests) 

We  are  interested  in  the  improvement  (reduction)  in  the  uncertainty 
(Cj)  which  can  be  written  as  follows: 


Oj,  as  shown,  can  now  be  substituted  into  Equation  (3-1)  to  provide  the  econo- 
mic value  of  the  Scimple  information.  Since  a-  = <7g/Zn  , this  value  will  change 
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with  the  number  of  tests.  The  economic  benefit  from  the  tests  will  be  the 
economic  value  of  the  sample  (test)  information  minus  the  cost  of  the  tests. 

For  example,  assume  the  following 

= 1 (reference  level  for  the  design) 

0^  = .056  (estimated  by  evaluating  the  uncertainty  range  in  the  mean 
estimate  after  reviewing  the  results  of  MIGHTY  EPIC) 

^ _ 5j  00_0, 000  million  material  cost/100%  change  in  strength) 

= .10  (typical  from  fragility  curves  in  Reference  1). 


From  Equation  (1),  EOL  = $111,704.  To  find  the  expected  average  gain  from 
the  test  information,  use  the  following  formulation. 


E [Value  of  Test  Information]  = EVTI 

(x-u^) 


dC 

dx 


/2ff  °1 


(x-u^) 


e 


dx 


dx  I 


dc  Jo  / °o^ 


(3-4) 


If  the  cost  per  additional  test  (cylinder)  is  = $50,000  then  the 
net  gain  from  additional  testing  is 


Net  gain  from 
additional  tests 


EVTI (n) 


- n C 


T 


(3-5) 


On  substituting  the  assumed  values  of  etc.  into  equations 

(3-4)  and  (3-5),  it  is  possible  to  determine  the  net  gain  financially  from 
the  addition  of  structures  to  the  test.  The  results  for  this  example  are 
shown  in  Table  3-2 . It  can  be  seen  that  one  additional  structure  would 
give  a small  net  gain  over  no  structures  or  over  two  structures  which  would 
cost  too  much  compared  to  the  benefit.  By  changing  the  material  costs  of 
the  future  facilities,  it  is  possible  to  do  an  analysis  of  the  optimum 
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number  of  additional  structures  as  a function  of  the  cost  of  material  in  the 
future  structures.  This  is  shown  in  Figure  3-2.  Figure  3-3  shows  the  average 
net  financial  gain  due  to  using  optimum  number  of  additional  tests  as  a 
function  of  the  costs  of  materials  in  future  systems.  Obviously,  there  is  a 
considerable  gain  as  the  future  systems  become  more  and  more  expensive. 

From  this  brief  analysis,  it  appears  that,  at  the  very  most,  one 
additional  structure  with  a predicted  failure  probability  of  .5  would  be 
desirable  in  the  DIABLO  HAWK  tests.  However,  it  appears  that  the  net  benefit 
from  the  additional  structure  will  be  small,  and  would  not  merit  the  delaying 
of  the  start  of  the  tests,  which  would  add  even  further  costs* 


n 

EVTI(n) 

n C.|. 

EVTI(n)-n  C.|. 

0 

0 

0 

0 

■■ 

54,579 

50,000 

4,579 

H 

69,350 

100,000 

-30,650 
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I Table  3-2.  Net  Gain  from  Additional  Structures 
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4, 


EFFECTS  OF  CORRELATION 


4 . 1 BACKGROUND 

As  indicated  previously  (Reference  1)  correlation  is  a significant 
but  frequently  omitted  aspect  of  a statistical  evaluation.  There  are  two 
prime  sources  of  correlation  in  the  lined  tunnel  problem: 

• The  free-field  stress  level,  if  high,  will  be  high  everywhere 
and,  if  low,  will  be  low  everywhere.  There  may  be  some  slight 
variations  in  the  value  from  point  to  point  (at  the  same  range) , 
but  the  heavy  trend  will  exist. 

• Strength  of  materials  such  as  the  reinforced  concrete  will  also 
tend  to  be  heavily  correlated  from  structure  to  structure 
because  of  the  commonality  of  material  source  and  manufacturing 
and  aging  processes. 

The  correlation  of  strength  of  structures  and  the  correlation  of  the  free- 
field  stress  have  no  relationship  to  each  other.  Thus,  in  a simulation 
procedure,  one  random  variable  is  generated  for  the  pressure  (and  appropri- 
ately scaled  for  all  the  tunnels),  and  another  (independent)  random  variable 
is  generated  for  the  material  strength  and  likewise  used  for  all  structures. 

The  easiest  problem  to  consider  is  the  case  where  the  structures  are 
fully-correlated  for  all  the  random  (structural  property)  variables.  That 
is,  if  the  concrete  is  especially  strong,  it  will  be  strong  for  all  the 
test  tunnels.  Similarly,  if  the  steel  is  slightly  too  thick,  it  will  be 
too  thick  for  all  the  structures,  and  so  on  for  '^ult' 

This  assumption  of  full  correlation  (in  all  the  random  variables)  will 
clearly  give  an  upper  bound  on  the  effects  of  correlation.  The  actual 
physical  problem  clearly  lies  somewhere  between  fully  correlated  (considered 
herein)  and  totally  uncorrelated  (considered  previously  in  Reference  1 ) . 

The  generation  of  failure  probabilities  for  the  fully-correlated 
problem  involved  developing  a small  F0RTRAN  program,  which  is  outlined  in 
Figure  4-1.  Beginning  on  the  left  of  Figure  4-1,  one  has  two  independent 
random  number  generators.  Since  (on  each  fragility  curve)  lies  between 
0 and  1,  the  procedure  is  to  generate  a random  number  between  [0,1).  The 
random  number  has  a uniform  probability  distribution  on  the  [0,1]  interval. 
Using  the  random  number  (called  "RANDY"  in  the  F0RTRAN  program)  one  enters 
the  fragility  curve*  and  reads  off  a structural  capability,  call  it  (o/o^) . 

_ 

Actually,  there  are  four  slightly  different  fragility  curves  in  the  program, 
corresponding  to  different  concrete  thicknesses.  See  Figures  4-5  through 
4-8  of  Reference  1 . 
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similarly,  one  generates  a log-normal  pressure  distribution,  cor- 
responding to  a dimensionless  range  R/R^  = 1.0.  Then,  knowing  the  range  of 

the  structure , one  proceeds  to  scale  the  pressure  variable  to  give  (P/P  ) 

th  ^ 

at  the  I—  structure.  Not  shown  in  the  Figure  is  a DO-LOOP  (over  I,  where 

I ranges  from  1 to  the  total  number  of  structures,  NSTRUC) . Within  the 

DO-LOOP  a test  is  made  to  see  if 

(p/pj^  > (4-1) 

If  equation  (4-1)  is  true,  then  the  I—  structure  will  fail,  and  a counter 
is  incremented.  For  example,  in  a single  pass  through  the  Monte  Carlo 
Loop,  5 structures  might  " fail " and  9 survive.  In  this  case,  the  variable 
FAIL  (5)  is  incremented  by  one. 

The  result  of  this ^procedure  (when  repeated  many  times)  is  to  popu- 
late a frequency  distribution  such  as  shown  in  Figure  4-2.  Note  that  on 
each  pass  through  the  Monte  Carlo  Loop,  the  random  number  generated  for  the 

structure  (P,)  is  used  for  all  the  tunnels  (hence  the  tunnels  are  100  per 

^ -2 
cent  correlated)  and  the  random  pressure  is  scaled  (using  P ~ R ) for  all 

ranges  (hence  the  pressure  loading  is  100  per  cent  correlated) . 

The  computer  program  was  relatively  straightforward,  and  it  is 
thought  to  be  error-free.  As  a check  problem,  one  range  [ (R/R^)  say]  was 
given  its  normal  value,  and  all  the  other  ranges  were  set  to  1000.  This 
check-out  procedure  corresponds  to  a single  structure  within  the  range  of 
the  weapon,  and  all  the  other  structures  out-of-range.  The  result  is  to 
produce  the  (uncorrelated)  failure  probability  of  a single  (independent) 
tunnel,  at  range  (R/R^)^^.  This  value  for  P^  was  then  checked  against  the 
curve  in  Figure  2-2,  and  agreement  was  obtained. 

After  thus  verifying  that  the  program  was  working  correctly,  several 
runs  were  made,  using  three  different  "factors  of  uncertainty"  in  the  log- 
normal pressure  distribution.  The  primary  result  of  correlating  the  struc- 
tures and  correlating  the  loads  is  that  the  probability  distributions 
(histograms)  begin  to  develop  significant  "tails"  and  the  standard  deviation 
increases.  This  result  is  in  agreement  with  physical  intuition  about  the 
problem,  since  "If  the  structures  are  all  extra-strong,  they'll  all  be  more 
likely  to  survive,  and  if  they  are  all  extra-weak,  they'll  all  be  more 
likely  to  fail."  Similar  arguments  for  emphasing  the  extremes  of  the 
probability  distributions  apply  when  the  correlation  of  the  pressure  is 
considered . 
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Individual  results  for  MIGHTY  EPIC  and  DIABLO  HAWK  are  discussed  in 
the  sections  which  follow. 

4.2  EFFECT  OF  CORRELATION  ON  DIABLO  HAWK  PREDICTIONS 

The  results  presented  in  Section  2.2  (cf.  Figure  2-4a  through  2-4c)  are 
based  upon  un-corr elated  random  variables.  Thus  the  concrete  strength  of 
one  cylinder  might  be  high  when  the  concrete  strength  of  its  neighbor  is 
low,  etc.  Similarly,  the  results  of  Section  2.2  assume  that  the  pressure 
loading  (i.e.,  free-field  stress)  is  also  random  from  cylinder-to-cylinder . 

From  a practical,  physical  standpoint,  the  assumption  that  the  cylin- 
ders are  all  uncorrelated  (with  regard  to  concrete  strength,  rock  properties, 
etc.)  is  not  correct.  Presumably,  the  same  construction  crew  made  each  of 
the  cylinders,  some  were  made  on  the  same  day,  they  were  made  using  the 
same  techniques,  etc.  Thus,  it  is  anticipated  that  certain  properties  of 
the  structures  will  be  highly  correlated.  Similarly,  the  pressure  loading 
(from  the  nuclear  event)  is  also  expected  to  be  highly  correlated,  such 
that  a high  free-field  st^ress  on  one  structure  implies  a high  free-field 
stress  on  its  neighbor,  etc. 

Note  that  the  pressure  loading  is  uncorrelated  with  the  material 
properties,  however.  That  is,  the  high  concrete  strengths  made  by  the 
construction  crew  are  totally  uncorrelated  with  the  effects  produced  by 
the  weapon.  From  a computational  standpoint,  the  simplest  case  to  consider 
is  that  where 

(i)  All  the  material  properties,  rock  strengths,  concrete  strength, 
etc.  are  fully-correlated  from  structure-to-structure, 

(ii)  The  pressure  loading  is  fully-correlated  from  structure-to- 
structure,  and 

(iii)  There  is  no  correlation  between  material  strengths,  cylinder 
properties,  and  the  applied  loading. 

This  problem  was  solved  computationally  by  using  two  independent 
random  number  generators  via  the  computer  program  outlined  in  Figure  4-1. 

The  results  of  these  calculations  are  shown  in  Figure  4-2  through  4-4, 
which  were  each  computed  for  a different  factor  of  uncertainty  in  the  free- 
field  pressure.  Referring  to  Figure  4-2  and  comparing  it  with  Figure  2-4a 
one  observes  that  the  effect  of  correlation  has  been  to  increase  the  prob- 
ability of  exactly  3 failures  and  exactly  8 failures.  Previously,  the 
probability  of  exactly  6 failures  was  quite  high  (see  Figure  2-4a) . Thus, 
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Figure  4-3.  Fully-Correlated  Histogram  for  Diablo  Hawk 
(Factor  of  2 Uncertainty  in  Pressure) 
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the  effect  of  correlation  has  been  to  modify  the  extremes  of  the  distribu- 
tion and  increase  the  standard  deviation. 

Similar  results  are  noted  when  Figure  2-4(b)  is  compared  with  Figure  4-3. 
Again  the  effect  of  correlation  has  been  to  modify  the  extremes  of  the  dis- 
tribution. This  effect  if  most  pronounced  for  the  larger  uncertainties 
in  the  free-field  pressure  (compare  Figure  2-4(c)  with  Figure  4-4).  With 
respect  to  the  predictions  for  DIABLO  HAWK,  since  the  effect  of  correlation 
has  been  to  widen  or  spread  out  the  distribution,  the  range  on  the  predicted 
number  of  failures  is  increased.  Thus,  although  in  Section  2.2  it  was 
predicted  that  "either  5 or  6 failures  would  occur,*  within  90  per  cent 
confidence  limits,"  the  fully-correlated  results  of  Figure  4-2  indicate 
that  between  3 and  8 structures  will  fail  in  the  DIABLO  HAWK  test,  within 
90  per  cent  confidence  limits.  Similarly,  for  the  higher  uncertainty 
factors,  the  range  on  the  number  of  predicted  failures  increases.  These 
effects  of  correlation  and  of  uncertainty  in  the  free-field  pressure  were 
also  demonstrated  when  calculations  were  preformed  for  the  MIGHTY  EPIC 
cylinders.  The  fully-correlated  results  for  MIGHTY  EPIC  are  discussed  in 
the  following  section. 

4.3  EFFECT  OF  CORRELATION  ON  THE  MIGHTY  EPIC  RESULTS 

Referring  to  Figure  4-5  (ta)cen  from  Reference  1 ) one  sees  that  the 
probability  of  4,  5,  or  6 failures  was  relatively  high,  when  no  correlation 
was  used.  Conversely,  Figure  4-6  shows  that  the  probability  of  exactly  5 
failures  is  reduced,  and  the  probability  of  exactly  8 or  10  failures  is 
increased,  when  the  structures  and  loading  are  fully-correlated. 

Previously,  without  correlations,  it  was  predicted  that  "from  4 to  6 
structures  will  fail,  within  90  per  cent  confidence  limits."  Again,  it 
is  noted  from  Figure  4-6  that  the  range  of  the  probable  number  of  cylinders 
failing  is  increased  when  correlation  is  used.  Thus,  from  Figure  4-6  one 
can  determine  that  the  probable  number  of  failures  lies  between  4 and  10, 
within  90  per  cent  confidence  limits,  for  a fully-correlated  MIGHTY  EPIC  test. 

With  higher  uncertainty  factors  (for  the  free-field  pressure)  the 
effect  of  correlation  is  to  emphasize  the  extremes  of  the  distributions, 
as  noted  previously.  (See  Figures  4-7  and  4-8,  for  example.) 

* 

For  a factor  of  1.1  uncertainty  in  pressure. 
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NUMBER  OF  STRUCTURES  WHICH  FAIL 

Figure  4-5.  Un-correlated  Histogram  for  Might  Epic 
(Factor  of  1.1  Uncertainty  in  Pressure) 
[Reference  1] 
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Figure  4-6 


Fully-Correlated  Histogram  for  Mighty  Epic 
(Factor  of  1.1  Uncertainty  in  Pressure) 
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4.4 


DISCUSSION  OF  CORRELATION 


It  can  be  argued  that  the  effects  of  correlation  should  be  such  as 
to  emphasise  the  extremes  of  the  probability  distributions,  as  was  dis- 
cussed in  Section  4.1.  In  brief,  if  the  properties  of  the  structures  are 
all  correlated,  then  they  are  more  likely  to  all  be  over-strength  (or  all 
under-strength)  and  all  survive  the  loading  (or  all  fail) . Similarly,  if 
the  pressure  loading  is  correlated,  then  it  is  more  likely  that  they  will 
all  survive  (if  the  pressure  is  low)  or  they  will  all  fail  (if  the  pressure 
is  high) . Thus,  from  a physical  standpoint,  it  is  understandable  why 
correlation  emphasizes  the  extremes  of  the  distribuitons . 

The  real  situation  in  the  (MIGHTY  EPIC/DIAGLO  HAWK)  structures 
experiment  falls  somewhere  between  the  "fully-correlated"  case  (treated 
herein)  and  the  "fully-uncorrelated"  case.  For  example,  the  cylindrical 
structures  all  were  made  by  the  same  contractor,  which  tends  to  make  them 
highly  correlated.  However,  the  rock  in  which  they  were  placed  may  vary 
from  one  location  to  another,  which  tends  to  make  the  properties  a little 
less  correlated. 

Sensitivity  studies  (Reference  5 ) have  shown  that  the  rock  proper- 
ties (like  the  ultimate  strength,  0^^^'  fracture  parameter,  k^^.) 

have  a significant  effect  on  the  structural  capability  (o/o^)  of  the 
tunnels.  Since  the  tunnels  are  fairly  widely  separated  in  the  rock,  one 
can  argue  that  these  sensitive  properties  are  relatively 

un-correlated,  and  hence  the  tunnel  strengths  should  be  uncorrelated. 
Conversely,  another  sensitive  parameter  is  the  concrete  strength,  f^'  , which 
is  expected  to  be  highly  correlated  (due  to  commonality  of  the  manufacturing 
process) . Thus,  there  cire  possible  arguments  for  high  correlation  among 
the  test  structures.  It  seems  clear  at  this  point  that  the  actual  struc- 
tures in  (MIGHTY  EPIC/DIABLO  HAWK)  fall  somewhere  between  the  two  extremes 
of  totally  uncorrelated  and  fully  correlated. 
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PROBABILITY  CALCULATIONS 
FOR  SPHERES 


In  addition  to  the  lined  tunnels  in  rock  (which  were  designed  by 
Merritt  CASES) , the  (MIGHTY  EPIC/DIABLO  HAWK)  tests  include  hollow  spherical 
concrete  structures  (designed  by  Agbabian  Associates) . Reference  6 
discusses  the  rationale  behind  the  test  design  for  the  concrete  spheres, 
and  Reference  7 presents  the  equations  used  for  designing  the  individual 
spheres  themselves.  It  is  noteworthy  that  several  fundamental  differences 
' exist  between  the  cylinder  design  (which  was  developed  by  Newmark,  Refer- 

ence 8 ) and  the  sphere  design  (which  is  based  on  Haynes's  work.  Reference  9). 

S 

I For  example,  Newmark 's  tunnel  design  procedure  involves  using  a 

I steel  liner  and  allows  plastic  deformation  of  the  steel,  the  surrounding 

? concrete,  and  the  grout/rock  medium.  Conversely,  the  sphere  design  uses 

no  steel  liner,  involves  fiber-reinforced  concrete  and  assumes  elastic 
! behavior  of  the  concrete  and  rock.  The  tunnel  design  by  Merritt  CASES 

i (Reference  5 ) accounted  for  uncertainties  by  using  relatively  standard 

1 design  practice  and  considered  maximum  and  minimum  values  of  rock  properties, 

concrete  properties,  etc.  On  the  other  hand,  Agbabian  Associates  used  a 
probabilistic  design  procedure  (based  upon  a paper  by  Ang  and  Cornell, 
Reference  10)  to  account  for  uncertainties  in  material  properties,  etc. 

Thus,  it  is  clear  that  many  basic  differences  exist  between  the  cylinder 
and  sphere  design  procedures.  These  differences  are  worth  examining  to 
increase  our  understanding  of  the  design  of  deep-based  structures. 

Other  questions  arise  with  respect  to  the  sphere  design  which 
1 concern  probabilistic  design  in  general.  The  results  presented  herein  and 

J in  Reference  1 are  based  on  a Monte  Carlo*  technique,  which  is  more 

[ general  and  can  treat  wider  variations  in  the  data  than  the  method  of  Ang 

[ and  Cornell  (Reference  10) . The  latter  approach  is  based  on  a perturbation 

procedure  which  contains  the  assumption  that 

« 1 (5-1) 

where  is  the  mean  value  of  the  i—  (uncertain)  parameter  (e.g..  Young's 
modulus  for  rock,  etc.) 

and  is  the  standard  deviation  of  the  i—  parameter. 

For  a discussion  of  the  term  "Monte  Carlo"  technique,  see  Reference  11. 
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Rather  than  go  into  a theoretical  comparison  of  the  probabilistic 
methods  used  by  the  various  contractors,  it  is  probably  of  more  value  to 
describe  the  individual  approaches  and  show  their  strengths  and  shortcomings. 
Figure  5-1  shows  three  methodologies  which  have  been  used  to  evaluate  the 
effects  of  uncertainties  in  structures  on  structure  strength.  The  first 
method,  I,  (Reference  7),  uses  a series  of  factors  multiplied  times  the  nom- 
inal strength  to  express  the  variation  in  the  actual  strength.  The  assump- 
tion is  made  that  the  uncertainty  in  these  factors  is  small  compared  to  the 
mean.  The  assumption  is  also  made  that  there  is  no  statistical  correlation 
between  factors.  The  approach  does  handle  systematic  error  by  the  incorpor- 
ation of  an  ignorance  uncertainty.  The  approach  works  only  for  small  uncer- 
tainties and  if  large  uncertainties  are  introduced,  will  produce  very  biased 
results.  Consequently,  this  approach  should  be  used  with  care,  and  always 
with  very  small  variations.  The  methodology  is  not  applicable  to  problems 
with  large  non-linearities.  A second  method,  II,  used  by  Merchant  (Reference 
12)  uses  a perturbation  technique  wherein  partial  derivatives  of  strength  are 
taken  with  respect  to  each  of  the  system  parameters.  A linear  statistical 
model  is  then  constructed  wherein  the  variance  of  the  strength  becomes  a 
function  of  the  products  of  the  squares  of  the  partial  derivatives  and  the 
variances  of  the  parameters.  The  model  as  used  did  not  contain  correlation, 
but  mathematically  is  capable  of  handling  correlation  between  the  parameters. 
Again,  this  model  is  valid  for  small  uncertainties  and  very  small  non-linear- 
ities, but  will  fail  with  large  uncertainties  and  non-linearities. 

The  third  appraoch.  III,  is  basically  a Monte  Carlo  method,  wherein  the 
model  is  always  preserved,  and  a simulation  is  used  wherein  sampling  from  the 
uncertainties  of  the  parameters  produces  a distribution  of  strengths  of  the 
structure.  This  method  is  the  most  exact,  but  is  also  the  most  expensive 
when  the  computational  costs  of  running  the  simulation  are  high.  From  the 
standpoint  of  comparison,  this  third  model  is  usually  used  in  checking  other 
models  such  as  the  Method  II  linear  statistical  model.  For  very  small  pertur- 
bations, Methods  II  and  III  agree.  There  has  been  no  direct  comparison  between 
Methods  I and  III  in  this  problem  to  substantiate  the  final  validity  of  Method 
I. 

Although  the  Monte  Carlo  technique  can  be  applied  to  the  sphere  design, 
and  fragility  curves  can  be  derived  for  spheres,  no  move  was  made  to  develop 
an  attendant  computer  program.  The  sphere  problem  was  studied  sufficiently 
to  outline  the  flow  of  the  calculations,  however,  and  the  results  are  given 
in  Appendix  A.  Another  problem  which  was  briefly  investigated  concerns  the 
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development  of  elastic-plastic  design  of  spheres,  which  is  discussed  in  the 
paragraphs  which  follow. 

The  Newmark  procedure  (described  in  Reference  8 ) allows  for  elastic- 
plastic  deformation  of  the  concrete  liner  and  the  rock.  Newmark  assumed 
that  the  tunnel  was  very  long  (the  plane  strain  approximation)  and  axi- 
symmetric  (no  variation  in  the  circumferential  direction,  0) . Consequently, 
the  tunnel  problem  is  taken  to  be  one-dimensional  - i.e.,  variations  along 
a radius,  r,  only  are  allowed.  Similar  approximations  (e.g.  , that  for  a 
sphere  the  loading  is  spherically-symmetric)  allow  the  sphere  problem  to 
be  treated  as  one-dimensional  (variations  allowed  in  the  radial  direction 
only) . Thus  the  question  of  whether  or  not  one  can  develop  a "spherically- 
symmetric  Newmark  procedure"  for  spheres  in  rock  naturally  arose. 

This  problem,  of  a spherical  cavity  in  an  infinite  elasto-plastic 
medium,  was  recently  analyzed  by  Durban  and  Baruch  (Reference  13) . Ref- 
erence 13  is  based  on  Durban's  DSc  thesis,  and  the  approach  is  not  as 
simplified  as  Newmark's  procedure.  However,  Durban  presented  a fully  non- 
linear analysis  (both  material  and  geometric  nonlinearity)  which  will 
simplify  somewhat  if  it  is  limited  to  small  strains  (e.g..  Eg  = .01  to  .05). 
Regardless  of  such  simplification,  Durban  gives  solutions  in  the  form  of 
closed  integrals  which  can  be  readily  evaluated  numerically.  The  extension 
of  these  results  to  include  a spherical  liner  in  a spherical  cavity  appears 
fairly  straightforward,  providing  one  is  willing  to  assume  perfect  bonding 
(continuity  of  stress  and  strain)  between  the  concrete  liner  and  the  rock. 

Although  the  subject  of  a "spherical  Newmark"  procedure  has  been 
examined  just  briefly,  three  results  stand  out: 

• The  spherically  symmetric  problem  is  somewhat  more  complicated 
than  the  cylindrically  symmetric  tunnel. 

• The  solution  is  amenable  to  numerical  integration  and  is  expected 
to  be  fairly  simple  to  implement  on  the  computer. 

• It  is  not  obvious  that  the  added  complexity  of  allowing  elastic- 
plastic  behavior  would  result  in  an  improved  design  (relative 

to  the  semi-empirical  procedures  which  Haynes  has  developed 
for  spheres;  see  Reference  9 ) . 

Further  work  to  answer  questions  about  the  elasto-plastic  design  of  spheres 
in  rock  was  beyond  the  scope  of  this  report. 
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CONCLUDING  REMARKS 


A primary  conslusion  of  this  study  is  that  5 or  6 cylindrical 
structures  are  expected  to  fail  in  the  DIABLO  HAWK  test,  with  5 tunnels 
being  heavily  damaged.  The  effects  of  initial  ovaling  (i.e.  initial 
imperfections)  from  the  MIGHTY  EPIC  test  are  thought  to  be  insignificant, 
however  this  estimate  remains  to  be  proven.  To  adjust  the  DIABLO  HAWK 
test  and  improve  the  probability  of  obtaining  better  design  data,  it  is 
tentatively  recommended  that  more  new  cylinders  be  added  to  the  test. 

The  effects  of  correlation  of  the  structural  strengths  and  of  the 
applied  pressure  loading  were  examined,  and  it  was  found  that  correlation 
emphasizes  the  extremes  of  the  probability  distrubutions . That  is,  the 
probability  of  a few  failures  increases,  and  the  probability  of  many  failures 
increases.  The  result  is  that  the  confidence  bounds  of  the  expected  number 
of  failures  increases.  For  example,  without  correlation,  5 to  6 failures 
are  predicted.  For  DIABLO  HAWK,  whereas  with  full  correlation  from  3 to  8 
failures  are  predicted.*  The  actual  physical  problem  will  lie  somewhere 
between  these  two  extremes  of  totally  uncorrelated  and  fully  correlated. 

The  study  showed  it  is  a relatively  straightforward  problem  to 
develop  fragility  curves  for  spheres,  but  time  and  money  constraints 
would  not  allow  their  development  herein.  A "spherical  Newmark"  procedure, 
for  elasto-plastic  design  of  spheres  in  rock,  might  also  be  developed,  but 
it  is  not  expected  to  be  as  simple  as  the  cylindrical  tunnel  procedure. 


♦These  predictions  assume  a factor  of  1.1  uncertainty  in  the  free-field 
pressure. 
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APPENDIX  A 

FLOW-CHART  FOR  SPHERE  CALCULATIONS 

Agbabian's  calculations  (Reference  7)  were  reviewed  to  determine  how 
fragility  curves  for  spheres  could  be  produced.  A possible  means  of  generat- 
ing the  fragility  curves  is  shown  in  Figure  A-1. 

Beginning  on  the  left  of  Figure  A-1,  one  selects  a free-field  pressure, 
P.  Then,  by  generating  independent  random  variables  and  using  equations 
which  relate  these  variables,  one  can  calculate  a (random)  in-plane  cracking 
pressure,  Next,  one  tests  to  see  if 

P > P , (A-1) 

- pi 

If  equation  (A-1)  is  satisfied,  a counter  is  incremented,  counting  the  fact 
that  failure  has  occurred.  Then  the  failure  probability  is  calculated,  and 
the  entire  process  repeated  in  a Monte  Carlo  Loop. 

The  result  is  that  for  a particular  pressure,  P,  there  will  be 
generated  a corresponding  failure  probability,  P^.  One  can  thus  develop  a 
fragility  curve,  P^(P),  of  failure  probability  as  a function  of  free-field 
pressure.  It  might  be  worthwhile  to  simultaneously  use  another  failure 
criteria  (say  implosion)  in  addition  to  in-plane  cracking.  Then  two  fragility 
curves  (one  for  each  failure  criterion)  would  be  developed.  (See  Figure  A-2) . 

Correlation  among  the  random  properties  and  the  pressure  loading 
could  also  be  treated,  following  the  guidelines  given  in  the  body  of  the 
report. 
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FAILURE 

PROBABILITY 


COUNT  THE 
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MONTE-CARLO 
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FAIL  = FAIL  + 1 


INCREMENT 

COUNTER 


THE  INTERACTION  FACTOR.  F, . 77-1:91 

IS  A FUNCTION  OF  K,  . E^ 

AND  R/t . SEE  REFERENCE  7 

Figure  A-1.  Flow  Chart  for  Generating  Fragility  Curves  for  Spheres 
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